"Pudding mold" band drives large thermopower in Na-jCoC^ 
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In the present study, we pin down the origin of the coexistence of the large thermopower and the 
large conductivity in Na^CoC^. It is revealed that not just the density of states (DOS), the effective 
mass, nor the band width, but the peculiar shape of the ai g band referred to as the "pudding mold" 
type, which consists of a dispersive portion and a somewhat fiat portion, is playing an important 
role in this phenomenon. The present study provides a new guiding principle for designing good 
thermoelectric materials. 



For the past decade, a cobaltate Na^CoC^ has been 
one of the most highlighted materials in the field 
of condensed matter physics in that it exhibits large 
thermopower [l| and magnetismjl] in the non-hydrated 
sodium rich systems, and superconductivity in the hy- 
drated sodium poor ones.[3| In particular, the discovery 
of large thermopower in Na^CoC^ [H and the findings 
in cobaltates/cobaltitesjl, [B|, @, 0, HI and rhodates^, |7(| 
that followed have brought up an interesting possibility 
of finding good thermoelectric materials that have rela- 
tively large conductivity. 

Theoretically, it has been proposed that the three fold 
degeneracy of the t^g orbitals as well as the strong elec- 
tron correlation effects plays an important role in these 
transition metal oxides. [1TJ Importance of spin degener- 
acy has been pointed out from experiments under mag- 
netic field. [l2j As for the orbital degeneracy, the t 2g or- 
bitals are split into an a lg orbital and two e' g orbitals 
due to the crystal field, and early first principles cal- 
culation predicted six hole pockets from the e' g bands 
in addition to the a\ g Fermi surface. (l3j However, an- 
gle resolved photoemission spectroscopy (ARPES) stud- 
ies [HI, HH, O [HI have revealed that the e' g bands lie 
O(100)meV~ O(1000)K below the Fermi level. There- 
fore, it seems unlikely that multiorbital effects affect the 
thermopower at least for T ~ O(100)K, where a large 
thermopower of S ~ O(100)/iV/K is already observed. [8( 
On the other hand, conventional calculations of the 
thermopower using the LDA bands already give large 
values. [ 13. Il4j suggesting that the band structure, not 
the electron correlation, is the key. In fact, it has been 
pointed out in some studies that the narrowness of the t^g 



bands is essential. [13|, |14|, LU| Although the narrowness 



of the bands must indeed be a factor, this alone cannot be 
the origin of the large thermopower because otherwise we 
would expect more good thermoelectric materials. Also, 
the problem of how the large thermopower and the large 
conductivity can coexist has not been addressed clearly 
so far. 

In the present study, we generally propose that a pe- 
culiar shape of the band (Fig[T]^a)), like the a± g band 
of the cobaltates (Fig[ljc)), can give large thermopower 



and large conductivity at the same time. In fact, in the 
preceding studies, we have proposed that the supercon- 
ductivity and the magnetism in Na^CoOo both originate 
from the peculiar shape of the a\ g band. [2(1 121| so adding 
the present study reveals that all of these interesting phe- 
nomena share the same origin. 

Our idea is the following. Using the Boltzmann's equa- 
tion approach, the thermopower is given as[22| 



(1) 



where e(< 0) is the electron charge, T is the temperature, 
tensors Kq and Ki are given by 



K„ = J2r(k)v(k)v(k) 
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Here, e(k) is the band dispersion, v(k) = V^e(fc) is the 
group velocity, r(fc) is the quasiparticle lifetime, /(e) is 
the Fermi distribution function, and /i is the chemical 
potential. Hereafter, we simply refer to (K. n ) xx as K n , 
and S xx — (l/eT)(Ki/K ) (for diagonal K ) as S. Using 
Kq, conductivity can be given as a xx — e 2 K^ = a = 1/p. 
Roughly speaking for a constant t, 



K ~ £>i 



,J B)i 



Ki ~ (kBT)Y±'(v% - v\) (3) 



(apart from a constant factor) stand, where S' is a sum- 
mation over the states in the range of \e(k) — fi\ <^ 
ksT, and va and vb are typical velocities for the states 
above and below [i, respectively. In usual metals, where 
va ~ vb, the positive and the negative contributions in 
K\ nearly cancel out to result in a small S (FigOJb)). 
Now, let us consider a band that has a somewhat flat 
portion at the top (or the bottom), which sharply bends 
into a highly dispersive portion below (above). We will 
refer to this band structure as the "pudding mold" type 
(FigOJa)). For this type of band with \i sitting near 
the bending point, v\ 3> v B holds for high enough tem- 
perature, so that the cancellation in K\ is less effective, 
resulting in \K\\ ~ (k B T)Y<'v 2 A and K ~ Y,'v\, and thus 
large \S\ - 0{k B /\e\) - O(100)/iV/K. A similar situa- 
tion where a large S originates from v\ 3> v B can be 
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FIG. 1: (color online) Schematic figure for (a) the pudding 
mold band and (b) a usual metal. The inset of (a) shows the 
top of the band shown in (c). The band dispersion for (c) 
t 2 /ti = -0.3, U/ti = -0.11 and (d) t 2 = t 3 = 0. The dashed 
line shows the position of fi at x = 0.7(n = 1.7) and T — 
300K. A contour plot of [v x {k)] 2 for x = 0.7 and T = 300K 
with (e) t 2 /ti = -0.3, ts/ti = -0.11 and (f) t 2 = ts = 0. 
Also region A (red outer half ring) and region B (green inner 
half circle or ring) are superposed. Only the k x > part of 
regions A and B is shown for clarity. 



realized for a simple parabolic band if [i lies very close 
to the band edge, but in that case, although the rela- 
tive ratio | is large, Kq is small because va, vb and 
the Fermi surface are small, so that the conductivity is 
small. By contrast, for the pudding mold band, the large 
va and the large Fermi surface (in 2D and 3D) result not 
only in large |^| oc S but also in large Kq oc a as well, 
being able to give a large power factor S 2 /p, which is 
important for device applications. (23| 

We now apply the idea to Na^OoC^. We take a 2D sin- 
gle band tight binding model H — — J2ij ^ij( c l c j + c \ c i) 
that considers hopping integrals between first, second 
and third nearest neighbor sites, t\, and £3. tijt\ — 
—0.3, t 3 /ti = —0.11 are introduced so as to reproduce 
the shape of the a\ g portion of the bands in the first 
principles calculation. [l3l [24j We take t\ = — 1000K 



as n =number of electrons/sites, and it is related to the 
sodium content x as n — 1 = x (as far as the e' g bands are 
fully filled). In FigQJc), we show the band dispersion, 
whose top around the T point (inset of FigQJa)) indeed 
has a form of a pudding mold. 

In the calculation of S, we first neglect the k depen- 
dence of t (so that r cancels out in Ki/K ). In Figj2£a), 
we show the temperature dependence of S for various 
values of band fillings. If we first focus on x — 0.7, we 
find that the overall temperature dependence observed 
experimentally for x ~ 0.71 in ref 8 as well as its large 
value of ~ lOO^iV/K at T ~ 300K is reproduced. In 
FigEJc), we show the Na content dependence of S at 
T = 300K, together with the experimental data of S 
at 300K taken from refH. Here again, we find a good 
agreement between the calculation and the experimental 
results. For comparison, we show in FigsOJb) and (c) 
the calculation result for t?. = £3 = 0, where the band 
top is just parabolic (FigOJd)). We see that the ther- 
mopower is strongly suppressed compared to the case 
with the pudding mold band. Note that for t% = t 3 = 0, 
the band width itself is nearly the same with the case of 
£2/^1 = —0.30 and t^/ti = —0.11, which indicates that 
the shape of the band is important. 

In order to see that our idea is working, we show 
in the bottom of FigQ] [ur(fc)] 2 , along with the regions 




to reproduce the band width obtained in the ARPES 
data, [II, Oil which about 60% of the LDA result 0, 0] 
due to electron correlation. The band filling n is defined 



FIG. 2: (color online) (a) Thermopower for t^/ti = —0.3 
and ta/ti = —0.11 plotted as functions of temperature for 
various x (denoted at the right edge), k dependence of r 
is not considered here, (b) Similar to (a) except t 2 — £3 = 
0. (c) Thermopower and (d) the normalized (normalized by 
the values for x — 0.7, t 2 /ti = —0.3 and tz/t\ = —0.11 ) 
resistivity (dashed) and power factor (solid) at T — 300K 
plotted as functions of x for t 2 /ti = —0.3 and tz/ti = —0.11 
(red) and t 2 — ts — (green). In (c), the experimental data 
at 300K (blue dots) are taken from ref.|J. 
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-2k B T < e(k) - p, < (region A) and < e(k) — fi < 
2k B T (region B). It can be seen that for the pudding mold 
band in (e), [v x {k)] 2 B /[v x {k)] 2 A <C 1, while for a parabolic 
band in (f) (t 2 = t 3 = 0), [v x (k)] 2 B /[v x (k)] 2 A ~ 1. It 
is not appropriate to say that the large S is due to the 
large DOS that originates from the flatness of the band 
top, because it is the dispersive portion (region A) that 
is positively contributing to S. 

The important expectation for the pudding mold band 
is not just the large thermopower, but also a relatively 
large conductivity and thus a large power factor. In 
FigEfd), we show the normalized resistivity pN and the 
power factor [S 2 /p]n with and without t 2 and £3 (here 
the x dependence of r is neglected) . The resistivity of the 
pudding mold band is only about a factor of ~ 2 larger 
than that for t 2 = £3 = 0, and this combined with the 
large thermopower results in a strong enhancement (10 
times larger at x = 0.7) of the power factor. 25 1 

So far, we have not focused on the effect of multiple 
Fermi surfaces (Fig[3la)), which has been proposed to 
be the ori gin of the superconductivity [20L [2H ] and the 
magnetism |2ll. [271 . l28j |. To see this effect, we plot in 
FigEIc) the t 3 dependence of S and [S 2 /p]n- Increas- 
ing £3 (or t 2 ) makes the local minimum structure at the 
r point deeper as shown in Fig{3fb) . S takes its maxi- 
mum at a certain £3, but the power factor continues to 
grow. This is because there are inner and outer Fermi 
surfaces both contributing to the conductivity, whose in- 
crease overcomes the decrease of S. If we go back to 
FigfS^d), this multiple Fermi surface effect is contributing 
to the low resistivity and high power factor for x ~ 0.8. 

In this context, it is interesting to estimate how large 
a power factor the multiple e' g hole pockets would give if 
they were present. In fact, the top of the e' g bands around 
the K point is again of the pudding mold type, with some 
corrugation that results in the multiple pockets. If we 
use the effective single band dispersion that resembles 
the upper portion of the e' g bands [2^], we find that the 
power factor is about two times larger than in the a\ g 
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FIG. 3: (color online) (a) Schematic plot of the pudding mold 
band with corrugation, p sitting at the dashed line results in 
inner and outer Fermi surfaces, (b) The variation of the band 
upon increasing ts/|ti| as — > 0.11 — > 0.24 with fc/ti = 
—0.3. (c) Thermopower (solid) and normalized power factor 
(dashed) plotted as functions of *3 / 1 ti | for t%/ti = —0.3, x = 
0.85 and T = 300K. 



case shown in FigfJ^d). Thus, if it were possible to push 
up the e' g bands so that the band top comes above the 
Fermi level for high Na content, the combined effect of 
both a\ g and e' g pudding mold bands may result in a very 
large power factor. Apart from this possibility, at high 
temperatures, and for relatively small Na content, where 
the Fermi level should sit close to the e' g band top, holes 
may be introduced at the top of the e' g band, which may 
be related to the enhancement of the thermopower above 
500K observed in ref.jl] 

We have seen that the consideration of the peculiar 
shape of the a\ g band alone suffices to understand the 
experimentally observed thermopower for x ~ 0.7 or T ~ 
300K. For x > 0.8 and T < 200K, however, this is not the 
whole story. In the calculated thermopower in Figf21[a), 
a hump structure around 100K found in the experiment 
for x > 0.8[H is not reproduced. In fact, this Na content 
and the temperature range is close to the region where 
the metallic magnetism, or the spin density wave (SDW), 
sets in (FigHKa)), whose origin has been shown to be 
due to the nesting between the inner and the outer por- 
tions of the Fermi surface. 21, 27| When the system is 



close to the SDW instability, the spin fluctuations become 
localized in q space around the nesting vector Q and in 
energy around to ~ 0. Thus, the quasiparticle scatterings 
by spin fluctuations occur with momentum transfer close 
to Q and with small energy transfer, so that, assuming 
a nearly isotropic Fermi surface, only the quasiparticles 
with \e(k) — p\ < Eb(T) is strongly scattered, where 
E B (T) is a certain energy scale that decreases with tem- 
perature (FigfJJb)). We leave microscopic treatment of 
this effect for future study, [13] but here, in order to qual- 
itatively take this effect into account, we assume a trial 
form of t as l/r(fc,T) oc A 2 {k, T)T+[1-A 2 (fc, T)]T 2 /^, 
A(k, T) = cosh -1 {[e(fc) - p{T)]/ E B {T)}. This form as- 
sumes that r in the strongly scattered regime is pro- 
portional to 1/T, while it crosses over to 1/T 2 in the 
weakly scattered regime. For E B (T), we assume a form 
E B {T) = k B T{T/&p + A Q (x), where £ 2 is the tempera- 
ture scale below which E B (T), apart from Ao, becomes 
smaller than the energy scale (~ k B T) relevant for Kq 
and K\ . We vary A with x as shown in FiglUJc) so as to 
reflect the tendency towards the SDW formation, which 
is the strongest around x — 0.8 ~ 0.85 experimentally. 
5, the Peltier conductivity a — Sa, and p (Fig(U(d)-(f)) 
calculated by using this r and taking 7 = 2, £1 = 1, and 
£2 = 0.1 (in units of |£i|) indeed reproduce the experi- 
mental results at least qualitatively. At low temperatures 
for x > 0.75, a exhibits a maximum and p crosses over 
from T to T 2 dependence in rough agreement with the 
experiment. Q This is because the quasiparticles some- 
what away from the Fermi surface have enhanced lifetime 
at low temperatures. [3~0| The thermopower S = a/a now 
has the hump structure at low temperatures for x > 0.8 
because K\ oc aT is enhanced more than Kq oc a since 
the former has larger contribution from the states away 
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terials that have pudding mold band(s) with the chem- 
ical potential lying close to its bending point. Roughly 
speaking, pudding mold bands tend to occur in geomet- 
rically frustrated lattice structures because there are sev- 
eral hopping contributions to s(k), which nearly cancel 
at some k, while not at other k. This, however, does 
not mean that "seemingly nonfrustrated" lattices are al- 
ways not good since, for example, a square lattice with 
sufficiently large diagonal hopping integrals (t') in addi- 
tion to vertical and horizontal hoppings (t) also results 
in a pudding mold band and a large thermopower (e.g., 
for t = -1 and t' = 0.45, S ~ 200^V/K at T = 0.3\t\ 
and n = 1.80). Another caution is that, even if a pudding 
mold band is present, a coexistence of a dispersive metal- 
lic band with a much larger Kg and a small K\ would re- 
suit in a suppression of Q2 bands K\™ d ) / (Z bands K b ™ d ) 
and thus a small total thermopower. Therefore, reliable 
band structure calculation is necessary in order to de- 
sign good thermoelectric materials based on the present 
guiding principle. 

Numerical calculations were performed at the facilities 
of the Supercomputer Center, ISSP, University of Tokyo. 
This study has been supported by Grants-in-Aid for Sci- 
entific Research from the Ministry of Education, Culture, 
Sports, Science and Technology of Japan, and from the 
Japan Society for the Promotion of Science. 



FIG. 4: (color online) (a) Schematic phase diagram concern- 
ing the relation between SDW and the thermopower. (b) The 
inset shows that when the spin fluctuations are localized in 
q, W space, the quasiparticle scatterings by spin fluctuations 
occur mainly with momentum transfer close to the nesting 
vector Q and with low energy transfer (solid arrow), while 
scattering with other q and large u (dashed arrow) is less ef- 
fective. Shown in the main panel is a schematic plot of the 
energy scale Eb(T) under which strong scatterings take place, 
(c) Variation of Ao with x. Calculation result for (d) the ther- 
mopower S, (e) the normalized Peltier conductivity a at and 
(f) the normalized resistivity pjv (normalized by the values at 
x = 0.7 and T = 300K) for various x using the trial form for 
r with 7 = 2, £i = 1 and £2 = 0.1 (in units of |ti|). 



from the Fermi surface due to the factor (e(k) — /i) in 
eq.([2]). At present, our form of r has no microscopic ba- 
sis, so we do not persist in these particular parameter 
values. In fact, for x — 0.85, e.g., by taking 7=1, 
£2 = 0.06, Ao = 0.08, a similar curve for S but with a 
smaller hump is obtained. 

To summarize, we have pinned down the origin of the 
coexistence of the large thermopower and the large con- 
ductivity in Na K Co02: the pudding mold a\ g band that 
has a dispersive portion and a flat portion with some cor- 
rugations at the top that can result in multiple Fermi sur- 
faces. The present study provides a new guiding principle 
for designing good thermoelectric materials: look for ma- 
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